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We present a construction of Hermitian operators and quantum states labelled by strings from a 
finite field. The distance between these operators or states is then simply related (typically, propor¬ 
tional) to the Hamming distance between their corresponding strings. This allows a straightforward 
application of classical coding theory to find arrangements of operators or states with a given dis¬ 
tance distribution. Using the simplex or extended Reed-Solomon code in our construction recovers 
the discrete Wigner function, which has important applications in quantum information theory. 


I. OVERVIEW 

Figure 1(a) depicts the binary Hamming cube - all binary strings of length 3 where strings that differ by one 
element are one edge length apart, strings differing by two elements are two edge lengths apart etc. The number of 
differing elements between two strings is the Hamming distance and finding useful arrangements of q -ary strings (with 
prescribed mutual Hamming distances) is the subject of classical coding theory. Let q = p n denote an integer that is 
a prime power. We will present a construction that associates q -ary strings with (i) Hermitian operators in Hilbert 
space of dimension dim("H) = q , and (ii) pure states in "H® 2 . We find a remarkably simple relationship between the 
Hamming distance of strings and the Hilbert-Schmidt or Fubini-Study distance of the corresponding operators or 
states, respectively. Because of the array of powerful coding-theoretic tools at our disposal, our construction may be 
useful for finding arrangements of quantum states or operators that would otherwise not be apparent. 



(a) (b) (c) 

FIG. 1. (a) Hamming cube for binary vectors of length 3, (b) the simplex code C C F| inscribed within it and (c) the facet 

operators A r corresponding to r £ Fo and their geometrical relationship to the Bloch ball (the subset of Hermitian operators 
corresponding to valid quantum states) for a single qubit. 

Figure 1(b) depicts the so-called simplex code, which is one member of family of q -ary codes that is well-defined for 
all prime powers q. Applying our construction to the codewords of this code (i.e., the vertices of the inscribed simplex) 
we find a set of operators that correspond to the phase point operators of Wootters’ discrete Wigner function. States 
that have non-negative quasi-probability representation in this Wigner function correspond to states that are not too 
far from any of the operators associated with the codewords of the code (see Fig. 2). This perspective on the Wigner 
function and its relationship with quantum state space may prove enlightening. 


II. MUBS AND FACE OPERATORS 

We will adopt the notation of quantum information theory so that the standard basis has elements | k) := £ 

C q = H and (-|-) is the inner product on H. Given an orthonormal basis B = {|0),..., \q — 1)}, a unit vector \v) is 
called unbiased if |(u|fc)| = W for all 0 < k < q — 1. We will focus on Hilbert spaces of prime power dimension q where 

it is known that q + 1 (the maximal possible number) mutually unbiased bases always exist. For non-prime-power 
dimensions the number of MUBs is lower bounded by the largest component in a prime decomposition of dim("H), 
but this is typically much lower than dirn('H) + 1. In subsequent sections we will be interested in connections between 
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MUBs and classical coding theory so we find it convenient to label MUB vectors with elements of the finite held, F g , 
containing dim('H) = q elements. In fact it is quite natural to use F g since many MUB constructions already use finite 
fields [1-3] so that e.g. Gauss sums can be used to prove the required overlap constraints. A complete set of MUBs 
has one more basis than the number of held elements so we label this basis with oo. 


Definition II. 1 Mutually unbiased bases: A complete set of MUBs in a Hilbert space of dimension dimfH ) = q is 
given by q +1 orthonormal bases {Boo, Bo, B±,...} = Use{oo f }&b, where each basis comprises Bb = {\^b),V G F g }, 
and overlaps obey 


I(V’bIV ! b')I 


— (1 — S b ,b’) + 8b,b'8v,v- 

Vd 


Cl) 


Our results do not depend on the specihcs of the mutually unbiased bases that we use. All that matters is their defining 
characteristic i.e., the pairwise inner products encapsulated in Eq. (1). In that sense, it is unnecessary that MUB 
vectors be labeled by elements of F g since any consistent labeling will do. We suggest that our choice is as convenient 
as any other and has the additional merit that the labeling is physically meaningful in at least one case, which we 
discuss in Sec. V. (In the context of quantum information this particular MUB construction is important because all 
basis vectors are eigenvectors of Pauli/Weyl-Heisenberg operators.) A good survey of different MUB constructions in 
power-of-prime dimensions is provided by Kantor [4]. Our construction also works without modification if we use a 
(necessarily incomplete) set of MUBs in dim("H) = p 2 that exclusively uses entangled basis vectors [5]. The association 
between Hermitian operators and F g -valued vectors is given by the following definitions, whose name derives from a 
geometrical interpretation described in Sec. V. 


Definition II.2 Facet operators: Using a complete set of mutually unbiased bases {\ifg),V G F q ,B G {oo,F g }} as 
in Definition II. 1, a facet operator indexed by a vector r G F 9+1 is defined as 

A r = l r ~' r = E WbWb B \~\- (2) 

B£F q ,oo 


For example in dim('H) = 3 a possible facet operator A r with r = [0,1, 2,0] corresponds to choosing the zeroth vector 
from the computational (Boo) basis, the first vector in the Bo basis, the second vector in the B\ basis and the zeroth 
vector in the B -2 basis. Dropping the requirement that we select a vector from every basis we arrive at the definition 
of a Face operator, 

Definition II.3 Face operators: Using a subset, of cardinality |r| (1 < |r| < q + 1), of a complete set of mutually 

\r\ 

unbiased bases, a face operator indexed by a vector r G F g , is defined as 

A r = £ |^?|-f |r| ~ Vg2 ~ g|r| + |r| ) V (3) 

BC{oo,F,} V 9 / 

|{B}|=M 


For example, we could drop the Boo and B\ bases from the previous example, and then A r with r = [ro,^] = [2,0] 
corresponds to taking the second vector from Bo and the zeroth vector from 62 . The definition for face operators 
completely subsumes the previous one since facet operators correspond to the special case |r| = q + 1. Nevertheless we 
have given them separate definitions as facet operators are the most interesting, and the simplification of the identity 
coefficient is not immediately apparent. 


III. FINITE FIELDS AND g-ARY CODES 

A field is a non-empty set F of elements with abelian addition and multiplication, satisfying the usual axioms e.g. 
distributivity. We denote as F g the finite held of order q = p n where p is a prime and n > 1 is an integer. The smallest 
number of times the unit element 1 gF 9 must be added to itself to produce 0 is the characteristic of the held, which 
is p, and consequently any element /3 G F g satishes p(3 = 0. If n = 1 and q = p then F g = Z p := {0,1,... ,p — 1} 
the integers modulo p. When n > 1 it is necessary to extend F p to F g with the addition of extra elements but we will 
not discuss the details of how this achieved. It will be sufficient to note that the nonzero elements of F g form a cyclic 
group of order q— 1 and a primitive element denoted a generates this whole group F g /{0} = {a, a 2 ,..., a 9-1 = 1}. 




3 


The 1-dimensional vector space F p n is also an n-dimensional vector space over F p . Let tr : ¥ q — p n h > F p be the trace 
map 


n— 1 

tr (f3):=J2P pk ( 4 ) 

fc=o 

then a standard result (useful in the context of Weyl Heisenberg operators later) is that for any 7 £ F 9 

w* r(/37) = qS 1} o where oj := exp(2ni/p). (5) 

/3eF„ 

A g-ary alphabet, that is a set of q distinct symbols, is naturally identified with elements of the finite field ¥ q . A 
word of length TV, w £ F^, is a string of TV symbols from F g and clearly there are q N distinct words of this fixed 
length. The most general definition of a g-ary code is as a subset C C F^ and the elements of C are called codewords 
(a good reference for all coding-related material is [6]). The Hamming distance 0 < A(v,w) < TV between two words 
v : w £ F^ is the number of positions in which v and w disagree. The Hamming distance is a metric on F^ so that 
expressions like A (it, w) < A(it, v) + A(v, w) hold. Using the Hamming distance we can define a ball/sphere of radius 
r around any word w via {u £ ¥^\A(v,w) < r}. Roughly speaking, good codes consist of codewords C C F^ where 
each codeword is the center of relatively large Hamming sphere, and this set of Hamming spheres fill the whole space 
without intersecting one another. The minimum distance d(C) of a code is given by d(C) = min{A(x, y)\x ^ y £ C}, 
and this is related to the radius of the empty Hamming spheres around each codeword. Two codes C and C are 
equivalent if they are related by trivial operations like permuting symbols or positions of codewords in a consistent 
way. Codes can be either linear or non-linear with the former typically being more amenable to analysis and simple 
encoding procedures. A linear code of length TV has q k codewords for some integer k > 0 and is denoted [TV, k, d\, 
whereas a nonlinear code has M codewords and is denoted (TV, M, d). From a purely combinatorial point of view, 
linear codes may be outperformed by nonlinear codes. 

The Hamming bound says that a g-ary code of block length TV and distance d has a cardinality |C| that is upper 
bounded by following expression 


r.rf-ti 

/TV\ 

Hamming Bound: |C| < q N / ^ . (g - 1)\ 

»=o 


( 6 ) 


and codes that saturate this bound are perfect e.g., the Hamming codes mentioned later. The Singleton bound says 
that a code C of block length TV and minimum distance d over a g-ary alphabet obeys 

Singleton Bound: \C\<q N ~ d+1 , (7) 

and codes that saturate this are maximum distance separable [7] (MDS) e.g., the simplex codes mentioned later. 

The standard notation for the number of codewords of Hamming weight i from the all zero codeword is 

Ai = |{tu £ F^|A(w,0) = i}|. (8) 

and it should not be confused with a face operator (the subscript and context should avoid this issue). The set 
{Ai|0 < i < TV} is the weight distribution of the code and is calculable using powerful tools like weight enumerators. 
Clearly for an (TV, M, d) code ^ = 

A code defines a vector space if and only if it is a linear code. A linear code encoding k units of information is 
described by a generator matrix G : F£ i —> F^ e.g., the Simplex code depicted in Fig. 1 has a generator matrix 


G 


simplex 


10 1 ' 


9i 

0 1 1 


32 


^simplex — ^9^1 Cl : b G F 2 }, 


= {( 0 , 0 , 0 ),( 1 , 0 , 1 ),( 0 , 1 , 1 ),( 1 , 1 , 0 )}. 


(9) 

( 10 ) 

( 11 ) 


The simplex code is well defined for all prime powers q and for all lengths of the form TV = (g m — 1 )/(g — 1) 
with parameters [TV = (g m — l)/(g — 1), fe = m,d = g 1 ”” 1 ]. The maximum length of a code that we may use in 
our construction corresponds to m = 2 and we will often refer to this code as the simplex code. This simplex code 
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sometimes goes by the name (doubly) extended Reed-Solomon code. In any event, our simplex code has generator 
matrix (recall that a is a primitive element of F g ) 


G 


simplex 


10 a a 2 ■■■ a g 1 


Si 

0 1111 1 


92 


^simplex — “1“ ^721^? b G Fg}. 


( 12 ) 

(13) 


The simplex code saturates the Singleton bound for all q but only saturates the Hamming bound for q = 3 where the 
simplex code is equivalent to the Hamming code 


G 


Hamming 


10 12 
0 111 


{q = 3). 


(14) 


The fact that this q = 3 code is doubly optimal (both MDS and perfect) arises from the following fact: the simplex 
code is dual to the Hamming code for all q but these codes coincide (the code is self-dual) for q = 3. The Hamming 
construction describes a family of codes with parameters [N = ( q m — l)/(g — 1), k = N — m, d = 3] so once again we 
consider m = 2 to describe the Hamming code for our purposes. This has a generator matrix with k = q — 1 rows i.e., 



' l 

0 

0 .. 

-ofl- 1 

-a 9 " 1 ' 



0 

1 

0 .. 

-a 9 " 1 

-a 9 " 2 


^Hamming 

0 

0 

1 .. 



(15) 





—a 9-1 

-a 2 






-a 9 " 1 

—a 



For any linear code C we can define an equivalent code C = C + w by adding a constant offset vector w to each 
codeword so that both codes have the same distance distribution. A standard coding technique, typically used for 
decoding, is to partition into cosets of a linear code, where each coset is identified (non-uniquely) by a coset leader 
w. This decomposition is depicted as a standard or Slepian array as in Table I where we have given an example using 
the binary simplex code of Fig. 1(b). 


Coset leader w 

Remainder of C + w 

(0,0,0) 

(1,0,1) 

(0,1,1) 

(1,1,0) 

(0,0,1) 

(1,0,0) 

(0,1,0) 

(1,1,1) 


TABLE I. Slepian array partitioning F| into cosets of the binary simplex code (11). The top row corresponds to the vertices of 
the tetrahedron in Fig. 1(b), whereas the second row consists of the same strings translated by (0, 0,1). Together the simplex 
code and its translate exhaust all 8 points of the binary Hamming cube. 


IV. DISTANCES IN QUANTUM STATE SPACE 

If we start with an operator of the form 


A r = Y,^B B w B B \-Ki q 


BeB 


then a fairly straightforward counting argument shows that 

Tr(A r ) = |r| - qK, 

Tr ((A r ) 2 ) = (g ~ 1 + M)|r| - 2|r|A' + qK 2 , 

= q when K = M ± ^ ~ gM + |r| 


(16) 


(17) 

(18) 

(19) 


where the last line explains the somewhat peculiar choice of identity coefficient that we adopted in Def II.3. Observe 
that face operators are clearly Hermitian A r = (A r )l since each term in the sum is manifestly so. We will examine 
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the geometrical relationship between these face operators and it is assumed that the same bases are used in the 
construction of two face operators A r and A 3 . The Hilbert-Schmidt inner product between these operators has 
remarkably simple expression, which is arguably the key insight of this work: 

Lemma IV. 1 Let A r and A s be face operators of the form (3), in a Hilbert space of dimension dimfH) = q , then 

Tr (A r A s ) = q — A(r, s) (20) 

where A(r, s) denotes the Hamming distance (number of differing elements) between vectors r, s £ Fg L 
Proof Insert the face operator definition from Eq. (3) and use the definition of mutually unbiased bases i.e., 

KV’bIV’b')! 2 = -(1 - $B,B') + Sb,B’Sv,V' (21) 

along with the fact that Yj, Sj = M — A(r, s). 1 

Since the operators A r are elements of the space of bounded linear operators, then the distance between two such 
operators can be characterized by the Hilbert-Schmidt metric. 

Corollary IV.2 The Hilbert-Schmidt distance between two face operators A r and A s of the form (3) is 

Dhs {A r , A s ) := ^Tr^-A*)^-^)] = ^2 [q — Tr (HM S )] = y/2A(r, s) (27) 

|r| • 

We can also identify normalized pure quantum states with vectors r £ FJj by using the Jamiolkowski isomorphism 
[8, 9] and the distance between quantum states is once again simply related to the Hamming distance, 

Corollary IV.3 Let A r and A s be face operators of the form (3), in a Hilbert space of dimension dim(TL) = q, then 
pure states \ J r ) £ C 9 given by \ J r ) = (I® A r ) Y^kew \hh)/y/q have trace distance and Fubini-Study distance 

Dt r (| J r ), | J s » := s/\ ~ \(J r \J s )\ 2 = - V2gA(r, s) - A *(r,s) , (28) 

q 

Dfs ( \J r >, I J s )) ■■= V2 - 2\(J r \J s )\ = y/2(l - |1 - A(r, s)/q\) , (29) 

where the latter simplifies to Dfs = \/2 A(r, s)/q whenever A (r,s) < q. 

Proof First note that, although face operators are not unitary in general, the Jamiolkowski isomorph obtained by 
applying A to one half of a maximally entangled state produces a valid normalized pure state (which is not generally 
maximally entangled). This can be seen using (J r \J r ) = Tr ((A r ) 2 )/q = 1 and similarly 

<■ J r \J S ) = Tr (A r A 3 )/q = 1 - A(r, s)/q. (30) 

The simplex code is equidistant with constant distance A = q between codewords so that {|J r ),r £ C S i mp i ex } forms a 
complete orthonormal basis in C 9 . 

It is interesting to consider how evenly the set of states {\J r ),r £ F 9+1 } is distributed with resepect to the Haar 
measure. Finite sets of states approximating the uniform Haar measure are well studied and go by the name of 
state t-designs [10, 11] (where t > 1 is an integer that quantifies how good the approximation is). The complete 
set of mutually unbiased bases described in Def. II. 1 comprises a state 2-design. Numerical calculations suggest that 
the set {|J r ),r £ F 9+1 } provides a poor approximation to a Haar-uniform distribution of pure states in C 9 . For 
instance, the purity of the reduced state p\ in a bipartite system quantifies how entangled the bipartite state is via 


1 If we want our face operators to have unit trace we can solve for a more general form 


so that 


In that case we find 


A r = j V IV'b B XV's B I - Klq 

(22) 

Bet3 

Tr(A r ) = J\r\ - qK = 1 

(23) 

Tr((W) 2 )= J2(9 - 1 + H)|r| 2M JK + q K>=q 

(24) 

J-J q + \ K- -WM(5 + 1). 

(25) 

V in i 

Tr (A r A s ) -q 9 + 1 A (r, s) 

(26) 


and the Hilbert-Schmidt and Fubini-Study distance measures can be derived from this. 
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1 < Tr {p\) < 1 where the lower bound is saturated for maximally entangled states. A result due to Lubkin [12] states 
that a Haar-uniform distribution of bipartite pure states has average subsystem purity (Tr(p 2 ))Haar = 2 q/(q 2 + 1), 
whereas we find 


q = 3 : 


( Tr (Pl))F| 


(|)9+(|) 72 59 

3 4 81’ 


(31) 


which suggests that entangled states may be under-represented in {|J r ),r £ F^ +1 }. 

The existence of a Hamming code (15) with parameters [g + 1, q — 1, 3] means that for all prime power dimensions 
there exists a set of facet operators of size |{A}| = g 9_1 wherein any two elements obey 


D HS (A r ,A s )>V6- 


(32) 


For q = 2 the Hamming code is simply C = {(0, 0,0), (1,1,1)} and the facet operators correspond to opposite corners of 
a cube in the space of Hermitian operators as in Fig. 2. Using codewords of the Hamming code then the corresponding 
set of states obtained via Corollary IV. 3 obey 


D FS (\J r ),\J s )) > 



(33) 


As well as knowing the minimum distance d = 3 there exist powerful tools (e.g. weight enumerators [13]) for calculating 
the complete weight distribution (8) of codes such as this. In this way we can enumerate the number of states | J a ) at 
any given (discrete) distance from a particular reference state | J r ). 


V. MUBS AND FACET OPERATORS USING THE WEYL-HEISENBERG GROUP 

The starting point for Weyl-Heisenberg operators in a Hilbert space of prime power dimension dirn('H) = q are the 
operators 


X (x)\k) = \k + x), Z(z)\k) = uj tT ( kzS> \k) x, z, k £ F g , w := exp(2ni/p), 


(34) 


which compose as 


X(x)Z(z)X(x')Z(z') = oj trx ' z X{x + x')Z{z + z'). (35) 

A Weyl-Heisenberg (generalized Pauli) operator, indexed by x,z £ F 9 , is a product of these A' and Z operators. 
From the composition law we observe that two Weyl-Heisenberg operators commute if and only if tr(xz' — x'z) = 0. 
The Weyl-Heisenberg operators generate a group that, modulo its center, has order q 2 . Consider a maximal abelian 
subgroup of this Weyl-Heisenberg group. Then any state that is a simultaneous eigenvector of all elements of this 
subgroup is a stabilizer state. Gross [14] showed that in a Hilbert space of dimension q = p n there are exactly 
p" TIL, (p' + 1) distinct stabilizer states. Our MUB constructions below are comprised of basis vectors that are 
stabilizer states. 

For odd prime powers, it turns out be convenient to impose a particular phase on the Weyl-Heisenberg operators 
so that they form the Weyl-Heisenberg group D of order |D| = g 2 , 

D ={D X}Z := w tr ^ X(x)Z(z)\x, z, £ F 9 } ^with ^ = 2 _1 /3 where 2 -1 £ F^j , (36) 

where individual group elements act as D X:Z \k) = w tr ^ +fcz |fc + x). Projectors onto rank -1 eigenstates of Weyl- 
Heisenberg operators (i.e., stabilizer states) can be constructed as [15, 16] 

I^XV’II = 1 E ^- kV) D k , k B, (37) 

q fceF, 


so that 


V,B£W q . 


( 38 ) 
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The set of states obtained by varying Eq. (37) over all B,V £ ¥ q , along with the computational basis Boo = 
{10), 11),...} is a complete set of mutually unbiased bases. One can check that the explicit form is given by 


W’b) 


— J2 u tr ^ Bk2 ~ vk) \ k), 

^ fee f. 


(39) 


and this 


is recognizable as the Ivanovic MUB construction [2, 3]. Using the composition law Eq. (35) we can deduce 


D^ V Ml\Dl, z = \€o +xs Ml +x \, (40) 

D*A*bM>b\DI* = ~ E 0J^- kv) D X!Z D k , kB Dl z , (41) 

q fceF, 

= ~Y1 u tl{ - k{V - z+xB)) D k ,kB, (42) 

q fceF, 

= U^- z+xB )«~ z+xB \- (43) 


Therefore the image of a facet operator under conjugation by a Weyl-Heisenberg operator is 
D A r D t _ /ir+a?[l,0,a,CK 2 ,...,a 9-1 ] — z[0,l,l,l,...,l] 

= A r+X9l ~ Z92 with 


' 1 

0 

a a 2 

• • • a«" 1 ' 


9i 

0 

1 

1 1 

1 1 


9 2 


(44) 

(45) 


where g i and gi are the generators of the simplex code (a similar expression was already pointed out in the prime¬ 
dimensional case in [17]). This is a very convenient way of understanding the orbit of facet operators under conjugation 
by the Weyl-Heisenberg group. It is also useful to have such a concise expression for the stabilizer states involved in 
the construction of a facet operator (for example, such a decomposition was used in [18] to construct a witness for 
quantum contextuality). 

For even-prime-power dimension, i.e., n qubits, it has been noted [1, 2, 19] that an Ivanovic-type MUB construction 
(39) over F g= 2 »* will not work without modification. Instead we must move to a slightly more general structure, the 
Galois ring GR(4, n), which has 4 n elements and its associated Teichmiiller set T = {0,1, £, £ 2 ,..., £ 2 ’^ 2 } with 2™ 
elements. Each element g £ GR{ 4, n) can be written g = a + 2b with a,b £ T and the trace map tr : GR{ 4, n) H>• Z 4 
is defined via 


n— 1 

tr (<7 = a + 2b) = ^ a 2 + 2 b 2 . 
k =0 


With these definitions we arrive at a MUB construction that appears formally very similar to the odd-prime-power 
case (39) 

I ^b) = - 7 = w 4 r(Sfc )+ 2 tr(vfe) | k) w 4 := exp{2m/4) = i (46) 


except now our labels are elements of T rather than F q . For the purpose of investigating geometrical relationships 
between face operators, the distinction between T-valued vectors and F 2 ™-valued vectors is irrelevant. From the form 
of the MUB vectors in Eq. (46) we can identify them as stabilizer states [20], just as we had in the odd prime power 
case. For the even q case we do not know of a similarly concise expression for the orbit of Weyl-Heisenberg operators 
acting on facet operators as we had in Eq. (44) although it should be possible. The Weyl-Heisenberg orbit of any A r 
with r £ T q+1 creates a simplex code e.g. for q = 4 we have 

D XiZ A[°’°’-- 0] Dl, = {A r | tr(r) e(00000), (01111), (02222), (03333), (10123), (11032), (12301), (13210), 

(20231), (21320), (22013), (23102), (30312), (31203), (32130), (33021)}, (47) 

where we are using coordinates tr(r) = (tr(r - 00 ), tr(n),...) £ Z\ rather than r £ T 5 . 


VI. THE DISCRETE WIGNER FUNCTION 

It is possible to represent finite-dimensional quantum states as probability distributions over a phase space of discrete 
points. However, to recover all the predictions of quantum mechanics we must allow the probability distribution to 
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take on negative values, in other words we must use quasi-probabilities. Such descriptions are referred to as quasi¬ 
probability representations, the most famous of which is the Wigner function. Wootters introduced a method of 
constructing discrete Wigner functions based on finite fields wherein vectors from a complete set of MUBs were put 
in one-to-one correspondence with the lines of the affine plane AG(2,F g ) [21], The connection with our work is that 
Wigner function of state p at some point in phase space is given by the expectation Tr(yfp), where A is one of the facet 
operators in Def. II.2 and which go by the name phase point operators in the context of Wigner functions. The tools 
and terminology established in previous sections allow for an interesting interpretation of the relationship between 
phase point operators with each other and with quantum state space. 

Wootters’ discrete Wigner function (DWF) requires a set of q 2 trace-orthogonal phase point operators, which 
corresponds to a set of facet operators {A r } with pairwise Hamming distance A(r, s) = q. We know that the 
codewords of Csimpiex satisfy this constraint, as do the codewords of every translate C s ; m piex + w for fixed w G F 9+1 . 
In this way we can obtain g 9_1 different DWF by partitioning F 9+1 into g 9_1 cosets of the simplex code via a Slepian 
array as in Table I. This partitioning is a coding-theoretic restatement of the concept of q 9-1 different “quantum 
nets” [21]. Hereafter, we will refer to a particular definition of DWF by its coset w, and it is understood that the 
MUBs used in constructing facet operators are those of Sec. V. 

The Wigner function of p at the point (x, z) 6 F ? x F, in phase space is denoted W Xt -(p), and is defined via 

W x Jp) = 1 Tr(A w+X9l ~ Z92 p) (fixed w <E F 9+1 ) (48) 

’ Q g 

The quantity Tr (A w+X9l ~ Z92 p) is the Hilbert-Schmidt inner product between the operator A and the density matrix 
p , and so demanding that W(p) > 0 is constraining p to be close to A in some sense. In fact the constraints 
W XiZ (p) > 0 Vx,z £ F 9 describe a simplex in H with q 2 bounding facets (hence the name for facet operators). In 
general our construction gives 

Simplex code in Hamming space <—> Simplex in Hilbert space. (49) 

Another geometrical object of interest is the single-particle (q = p ) stabilizer polytope defined as the convex hull of 
all p(p + 1) stabilizer MUB vectors | tpg). Cormick et al. [22] showed that a halfspace description of the single-qudit 
stabilizer polytope is given by 


Stabilizer polytope := {p| Tr {pA r ) > 0,Vr € F][ +1 }. (50) 

From the discussion in the previous paragraph we see that the stabilizer polytope is the intersection of all simplices 
associated with the simplex code and all its cosets (see Figure 2 for an illustration of the p = 2 case). 



(a) (b) (c) (d) 


FIG. 2. (a) Applying the facet operator construction to the simplex code C s i mp i ex and (b) to a coset C s i mp i ex + w (c) The 
intersection of these two simplices produces the stabilizer polytope (octahedron) (d) See also [23-25] for related geometrical 
discussions. 

We are interested in non-negatively represented pure states, that is states | <j>) such that W XtZ {\(j)){(j)\) > 0 Vx, z £ ¥ q . 
In Figure 2(a) we see that nonegatively represented pure qubit states in the DWF with w = (0,0, 0) are those that are 
both (i) on the surface of the Bloch sphere, and (ii) contained within the tetrahedron. Figure 2(b) illustrates the same 
idea for the DWF with w = (0,0,1). Cormick et al. [22] showed that the only pure states that are non-negatively 
represented for all q 9_1 Wigner functions (simultaneously) are the q(q + 1) stabilizer MUB states used in the 
DWF construction and this is illustrated in Figure 2(c,d). When we move to qutrit (or any odd prime dimension) 
state space the story changes slightly from the qubit case we have depicted. There are now g 9_1 = 9 different DWF 
and we are guaranteed by [22] that all q(q + 1) = 12 stabilizer states \ipg) are non-negatively represented no matter 
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which DWF we use. However, a result by Gross [14] says that a pure state is non-negatively represented in the DWF 
with w = 0 if and only if it is a stabilizer state. For qubits the set of non-negatively represented pure states is of finite 
measure, but for odd-prime qudits it is exactly the set of p(p + 1) stabilizer states. Moving on to multiple particles of 
odd prime dimension (q is an odd prime power) then it seems that Gross’ choice of DWF is the unique one obeying 
the discrete version of Hudson’s Theorem [14] : a pure state is non-negatively represented if and only if it is a stabilizer 
state. A priori we know that at least q{q + 1) stabilizer states will be positively represented but this only represents 
an exponentially small (in n) fraction of all p n n" =1 (p z -(-1) stabilizer states. Hence the discrete Hudson theorem is a 
geometrically remarkable fact, as well as having practical relevance for questions surrounding fault-tolerant quantum 
computing [26, 27], resources theories [28] and contextuality [18]. To see that w = 0 recovers Gross’ choice of DWF 
insert (39) into yp"=[o.o, -.,o] anc j simplify to obtain J2k&w q l&)( — &| he., the discrete parity operator (the parity operator 
also forms the starting point for the continuous Wigner function [29] ). This particular instance of Wootters’ discrete 
Wigner function is also singled-out by its highly symmetric properties [21, 36, 37]. 

We hope that our way of analyzing these Wigner simplices and their relationships with each other and with the 
set of quantum states will prove enlightening. In principle, we could construct a Wigner-like representation using 
the Alltop [30, 31] MUB vectors |^°Vb) = ^ S/ceF u) tT ( ak3+ 5 Bk2 ~ vk ')\k), which are equal to the Ivanovic MUB 
vectors for a = 0 and unitarily equivalent but highly non-stabilizer [32, 33] otherwise (a ^ 0). One could also apply 
our Wigner simplex construction to unitarily-inequivalent MUB vectors [4]. Note that Bengtsson and Ericsson [38], 
without restricting to stabilizer MUBs, have studied the equivalent of Wigner simplices and the stabilizer polytope 
(the complementarity polytope) and their relationship to quantum state space. Their construction of a simplex via 
mutually orthogonal Latin squares is isomorphic to our simplex code construction [7] in dimension q = p n . 


VII. SUMMARY 

We identified a construction relating the Hamming distance between g-ary strings to the Hilbert Schmidt distance 
between certain Hermitian operators and the Fubini-Study distance between certain states. Any q -ary classical code 
of length up to N < q + 1 is suitable and our hope is that the ability to use a vast array of coding-theoretic tools 
(distance distributions, automorphisms etc.) will prove useful in the quantum context. The types of operators and 
states that our construction provides are somewhat limited so it is unlikely that our results are directly applicable to 
outstanding open problems like the existence of symmetric informationally-complete positive operator-valued measures 
(SIC-POVMs). 

One topic for which our results are certainly relevant is the discrete Wigner function, where the expectation value of 
our operators are quasi-probabilities representing quantum states. The discrete Wigner function is both foundationally 
interesting as well as practically relevant for fault-tolerant quantum computing [26-28]. We showed how a famous 
family of maximum distance separable codes, the simplex codes, when applied via our construction, reproduce the 
Wigner simplex in Hilbert space. More generally, our results represent a convenient tool for working with, and a novel 
way of thinking about, Wootters’ Wigner function in arbitrary prime-power dimension. 
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